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The mystery of high temperature superconductivity in layered cuprates is tightly related to 
their common patter of the doping-dependent evolution from the antiferromagnetic 
insulator at zero doping to the overdoped metal. A number of experimental and theoretical 
studies indicate that the evolution is not smooth and a critical point separates the 
underdoped (UD) and overdoped (OD) regions. It is tempting to associate such critical 
point with the pseudogap formation for p<p*, with p* = 0. 19-^0.24 13 . No doubt that the 

proximity of the pseudogap and the superconductivity with two energy scales, T* and T c , 
is essential for high-T superconductivity 4 . On the other hand, the Hall effect 
measurements suggest that the critical point and the quantum phase transition (QPT) take 
place at optimal doping, p opt =0A6 5 . To resolve this controversy, here we performed 

precise density of states calculations and found that there are two QPTs and the 
corresponding critical concentrations associated with the change of the Fermi surface 
topology upon doping. 

Angle-resolved photo-emission spectroscopy (ARPES) reveals a change of the Fermi 
surface (FS) topology from the small hole pockets to the large hole FS near the optimal doping 6 ' 
1 . This provides a link between the QPT and changes of the FS. Here we will apply the general 
Lifshitz ideas on the QPT induced by the FS transformations. But first of all, we will discuss 
how these transformations are induced by doping. 

It is easy to obtain the large FS in cuprates by a single-electron approach like the local 
density approximation (LDA) or the tight-binding method. However, to get the small hole 

pockets around the (±tt/2, ±n/2\ points of the Brillouin zone requires going beyond the weak- 
coupling approximations and taking the strong electronic correlations into account. Such small 
pockets have been found in a doped antiferromagnetic (AFM) Mott insulator by exact 
diagonalization 9 and Quantum Monte Carlo 10 calculations for the finite clusters as well as by a 
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perturbation treatment of the infinite lattice " . According to these studies, after the long-range 
AFM order vanishes with increasing hole concentration n h = 1 + p (in La 2 _ x Sr x Cu0 4 , p - x ) a 

short-range AFM order still persists even at optimal doping 14 . The short-range magnetic order 
determines the self-energy and hole dispersion resulting in the small hole pockets around the 

(+7r/2, +tt/2) points in the UD cuprates, its fluctuation results in the pseudogap formation 15 ~ 17 . 

Due to the strong electronic correlations intrinsic for cuprates a theory of the electron dynamics 
has to fulfill a "no-double occupancy" constraint. This constraint is introduced explicitly in the 
mean- field theory of a d -type superconductivity within the RVB approach for the t — J model 
18 , and in the variational Monte-Carlo studies 19 . 

Contrary to the phenomenological approaches, like assuming that the second order QPT 
exists at p = p c 20 , we deal with a microscopically derived t — t' — f — J* model without free 

parameters. All parameters, including intraband t 1 } and interband t l l hoppings, have been 
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calculated within a hybrid LDA+GTB approach , which combines the ab initio LDA and the 
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generalized tight-binding (GTB) method for strongly correlated systems. Both short-range 
magnetic order and "no-double occupancy" constraint are incorporated in our approach by the 
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Hubbard X -operators perturbation theory, which was developed earlier for p - d ' and t—J 

24 ' 25 models. The method of irreducible Green functions used in this theory is similar to the 
Mori-type projection technique, with the zero-order Green function given by the well-known 
Hubbard I approximation. Beyond it there are spin fluctuations. To provide a description of 
them, the self-energy was calculated in the noncrossing approximation by neglecting vertex 
renormalization that is equivalent to the self-consistent Born approximation (SCBA) 23 . 
Resulting electron self-energy contains the space-time dependent spin correlation function. Note 
that at low temperatures T <10 K the spin dynamics is much slower than the electron one. A 

typical spin fluctuation time, 1CT 9 sec, is much larger than the electronic time 1(T 13 sec 26 ; that is 
why we can safely neglect the time dependence of the spin correlation function. The self-energy 

becomes static, Z(k,&>) — >Z(k), and is determined by static spin correlation function 
C 0n =(SqS~^ and kinetic correlation function K 0n =^(cL c »rr) between sites and n (see the 

methods section). These correlation functions and their Fourier transforms C q and K q represent 

the AFM short-range order and the valence bond order, respectively. A self-consistent 
calculation of these functions up to n = 9 (ninth coordination sphere) together with the chemical 
potential as functions of doping x reveals that all C 0n decrease with doping. While the nearest 

neighbor spin correlation function C m is finite for all studied x up to x = 0.40 with a kink at 
x = p* = 0.24, more distant spin correlations fall down to zero for x > p*. 

The calculated FS twice changes its topology with doping . Small hole pockets around 
(+^/2, +^r/2) points are present at small doping; then they increase in size and touch each other 

in the non-symmetric points k = +;r(l, + 0.4) at x cl = p opt =0.151. Above p opt , there are two 
Fermi surfaces around {jt,7t) with outer being a hole-like and inner being an electron-like. The 
electron Fermi surface collapsed at x c2 = p* = 0.246 , and at x > p * we have only one large hole 
surface around {n,7i) . Similar conclusion on the coexistence of hole and electron FS at some 
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intermediate doping have been also drawn recently ' . 

From the previous consideration it follows that the FS topological transitions in cuprates 
are induced by doping and they are due to the non-rigid band behaviour of the quasiparticles in 
the strongly correlated systems. According to the general Lifshitz analysis , for the three 
dimensional (3D) system a change of topology at the energy s = s c either by appearance of a 

new segment (like we found at p * ) or by change of its connectivity (like at p ) would result in 
the additional density of states, 8N{s) ~ (e-£ c ) 1 ' 2 , and the change in the thermodynamic 

potential, SD. ~ (s F -£ c ) 5 ' 2 (the QPT of the 2.5-order), where s F is the Fermi energy. However, 

due to the strong anisotropy of electronic and magnetic properties, cuprates are quasi-2D and not 
isotropic 3D systems. The electron hopping perpendicular to the Cu0 2 layers in a single-layer 

La 2 _ x Sr x Cu0 4 (LSCO), Bi 2 Sr 2 Cu0 6+s (Bi2201), etc. is negligibly small. We do not consider here 
YBa 2 Cu 3 7 _ s (YBCO) and Bi 2 Sr 2 CaCu 2 O s+s (Bi2212) with two Cu0 2 layers in the unit cell 
where the bilayer splitting of the Fermi surface appears, so we calculate the density of states 
(DOS) for the electrons in the doped single Cu0 2 layer. 

The change of the FS topology at x cl = p opt results in the logarithmic divergence of DOS 
(Fig. 1), while the emergence of the new electron-like pocket below x c2 = p* results in a step in 
DOS (Fig. 2). The total DOS is a sum of the singular and regular contributions. We would like to 
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stress that both log and step DOS singularities are in perfect agreement with the general 
properties of the van Hole singularities for the 2D electrons. Contrary to the 3D systems, the 

thermodynamical potential for the 2D electrons has a singular contribution <5n.~(e F -£ c f for 

the step singularity and <3T2 ~ {s F —£ c ) 2 \n\s F -s c \ for the log singularity 29 . Thus QPT at 
x c2 = p* is of the second order, while at x cl = p opt the singularity is stronger. It is immediately 
follows that the Sommerfeld parameter y in the electronic heat capacity y = CjT has also a 
singular step contribution at x < p* , and 8y cc \n(s F -£ c ) °c ln|jc-x o ^| near x cl = p opt . 

To check whether the coincidence of x cl with p and x c2 with p * is occasional or not 
we have calculated the superconducting critical temperature dependence T c (x) in the same 
model 30 and the kinetic energy as a function of doping. The T c (x) dependence is an inverse 
parabola with the maximum at x opt (see inset in Fig. 1), which indeed equal to x cl . Note that it is 
not a coincidence since like in the BCS theory the maximum in T c (x) is determined by the 
maximum density of states, and at x cl we have a log singularity. Kinetic energy, 

E ki„=zl t 0n K 0n> C 1 ) 

n 

reveals a remarkable kink at x c2 =p* (Fig. 3). Above p* , E kin (p)/E kin (p*)U l + p that is 

expected for a conventional 2D metal with hole concentration n h = 1 + p and E Un ~ s F ~ n h . The 

extrapolation of this law below p * (shown in Fig. 3 by the blue dashed line) reveals that actual 

E kin is smaller. We associate this depletion of the kinetic energy with the pseudogap formation 

and try to fit it with a simple free electron gas with a triangular pseudogap DOS (Loram-Cooper 
model ): 



N(e) = 
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Here E — J(p*—p)/p* is a doping dependent pseudogap, J is the exchange parameter. This 

fitting is shown by the filled triangles and it reflects some decrease of the kinetic energy due to 
the pseudogap but do not provide a quantitative agreement. Apparently, better fitting is given by 
the exponentional law, 

E ka {p)lE kin {P*) = ^p{-4E g (p)/j} . (3) 
Anyway, this analysis confirms that the QPT at x c2 is indeed related to the pseudogap and the 
coincidence of x c2 and p * is not occasional. 

A singular contribution to the Hall coefficient near the optimal doping has been measured 
for Bi 2 Sr 51 La 0A9 CuO 6+s single crystals and for La 2 _ x Sr x Cu0 4 thin films 5 under strong magnetic 

field of 60T. According to our theory, these extra carriers are induced by the singular DOS. To 
continue discussion in terms of the critical points, not the critical energies, we note that near the 

critical point e F (x)-£ cl = kyx-x opt j . In Fig. 4 we plot the calculated singular DOS A^ sing (z), 

z = x- x opt , together with the singular contribution to the Hall coefficient, 

n Haii (1 • )~ n Haii (10QK") ~ '. The optimal doping in the LSCO thin film x opt = 0.17 is shifted from 

the bulk value x opt =0.15 in the Bi2201 which may be due to the strains in the films. The 

general agreement of the calculated singular DOS and Hall data provides further support for our 
analysis. 
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In summary, we have shown that there are two critical points in the cuprate's doping 
dependence. The first one is related to the change of the FS connectivity and logarithmic 
divergences of DOS and of electronic heat capacity parameter y at the optimal doping 

Popt =0.151. Also, we associate this QPT with the experimentally observed singular doping 

dependence of the Hall coefficient 5 . Moreover, the logarithmic enhancement of DOS leads to 
the maximum in the doping dependence of superconducting critical temperature T c at the same 

critical point x = p t . This is in agreement with the previous calculations of T c for the magnetic 

mechanism of d 2 2 pairing. The second QPT is associated with the collapse of the electron-like 

x —y 

FS pocket at /? — > p* = 0.246 and results in the step singularities in DOS and in Sommerfeld 
parameter y . We have found the depletion of the hole's kinetic energy below p* and ascribe it 
to the pseudogap formation at p< p* . Thus the two energy scales in cuprates measured by T c 

and T* are both related to the QPTs and evolution of the cuprate's electron structure with 
doping. The very existence of both logarithm and step singularities in DOS are in perfect 
agreement with the general properties of the van Hole singularities for the 2D electron systems. 
But the concentrations of doping at which these singularities approach the Fermi level and start 
to govern the behavior of the system are determined by the strong electronic correlations and 
scattering on the associated short-range AFM order. 

"Methods" 

To properly fulfill the "no-double occupancy" constraint at every step of our calculations 
we use the Hubbard X -operators, X pq = \p)(q\, where \p) and \q) are the local eigenvectors 

corresponding to three states: one-hole state \<j) , <j = ±\/2, and the Zhang-Rice singlet, \S}, 
which is a two-hole (zero-electron) state. Constraint results in the exclusion of two-electron 
(zero-hole) state 1 0) = d 10 p 6 from the local Hilbert space. In the Hubbard model, there are two 

fermions with spin a : X 0a and X sS (a = -a) , forming the lower (LHB) and the upper (UHB) 

Hubbard bands. In the hole representation the LHB corresponds to an empty conduction band of 
electrons and the UHB corresponds to a valence band. The "no-double occupancy" constraint 
means the absence of the direct excitations from and to the LHB; the doped holes occupied only 
the UHB. Nevertheless the virtual interband (between the LHB and the UHB) hopping t n results 

in the exchange interaction J fg ={tf g ) /u efS . For La 2 _ x Sr x Cu0 4 , all intraband and interband 

hopping parameters (t" g and t 12 g ), single-site energies of holes in p - and d -orbitals, and the 

charge transfer gap U eff have been calculated by the ab initio LDA+GTB method 10 . The low 

energy effective model is the t — t'—f — J* model where J* means that besides the Heisenberg 
exchange term a three-site correlated hopping H 3 is also included. In terms of the Hubbard X - 

operators the Hamiltonian H t _ Jt = H t _j + H 3 is given by 

H t _ J =^{s-»)X 7 + ^2{s-v)Xf + X tlXfXf +1/4 x Jf g { X T X T ~ X T X 7 )>( 4 ) 

t 12 t 12 

j i \ fin mg ( v Sa v au v aS \rScr v aS V S2 \ 

H l = L Jj~\ X f X ™ X S ~ X f X m X g )■ 

Here hole creation operator is now c\ = 2a X Sa and its algebra is different from the bare 
fermion's one. The spin operators are also easily expressed via X -operators, S + = X aa , 
S z =i y X aa -X ffff )/2. The self-energy calculated in the SCBA 23 is determined by the spin 

correlation function C(q,a>) and results in the finite quasiparticle lifetime, ImZ(A:,ii) ^0. As 
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discussed above, at low temperatures we can neglect spin dynamics, C(q,co) — > C q , and then we 
have ImZ = 0. Nevertheless the doping evolution of the Fermi surface is determined by the 



ReZ , and it is qualitatively similar in our approach and in the approach which properly takes 



ImE into account 23 . The Green function for a hole moving against the background of short- 
range AFM order is written as 



G{KE) = ({x° 2 \x' k 
where 

V ' \ + xN^ 



l + x 



E-{s q +U^-m) — *, 
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respectively. In contrast to 



^eff 2 ^eff 
Here, t k and F k are the Fourier transforms of hoppings t 1 ! and t l i , 
approach of Ref. 23, we calculate correlation functions C q and K q self-consistently together 
with the chemical potential. To get the spin correlation function we also obtain the spin-spin 
Green function ((^f |^T)) 

in a spherically-symmetric spin liquid state 24 with (S 12 ) - and 

the equal correlation functions for each spin component, (SqS~\ = 2{slS z n \ = C 0n . Both C 0n and 
K 0n are essentially doping-dependent . 

It should be stressed that the standard DOS calculations with routine precision (400x400 
points in the quarter of BZ) we used before to solve the T c equation for magnetic mechanism of 

d x2 _ y2 -wave pairing is not enough to find the effect of QPT on DOS. To get the results 
presented here we used 10 5 xl0 5 k -points which lead to the increase of precision by 625 times. 
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Fig. 1. Regular, singular, and total density of states N(s—s cl ) near the optimal doping 

e c\ =£ f{ Popt ) as calculated from the Green function (5). Dotted line shows the log fitting. In 
the inset, the doping dependence of the superconducting critical temperature T c (x) is shown; the 

optimal doping is 0.151 . 
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Fig. 2. Regular (hole pocket), singular (electron pocket), and total density of states 
N(s-£ cl ) near the pseudogap critical point e c2 =s F (p*). Below p* = 0.24 (s < e c2 ) a 
singular step-like contribution to the total DOS appears. 
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Fig. 3. The doping dependence of the dimensionless kinetic energy E hn {p)/E Un {p*). 

Calculated dependence shown by the filled red circles. Above p * it obeys a conventional law 

and proportional to (l + p^ . The extrapolation of this law to the region p < p* (blue dashed 

line) emphasizes the depletion of part of the kinetic energy in the pseudogap region. Calculation 
for the idealized triangular pseudogap model (2) is shown by the filled triangles. 
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Fig. 4. Comparison of experimental singular Hall coefficient for Bi 2 Sr 05l La 049 CuO 6+s (a) 
and for La 2 ^ x Sr x Cu0 4 (b) and our calculated singular DOS N siDg [s F (x)) near the optimal 

doping. Agreement with results on bulk single crystals (a) is better than with results on thin 

films (b). 
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